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Abstract 

We give a brief overview of the results obtained in [TJ, concerning the rate 
of energy loss of an accelerating quark in strongly-coupled N = 4 super- Yang- 
Mills, both at zero and finite temperature. For phenomenological purposes, 
our main result is that, when a quark is created within the plasma together 
with its corresponding antiquark, the quark starts feeling the plasma only after 
the q-q separation becomes larger than the (w-dependent) screening length, and 
from this point on the motion is correctly described by the analytic energy loss 
formula previously derived by Herzog et al. and Gubser within the stationary 
or late-time approximations. The present text is a slightly expanded version of 
two talks given at the XIII Mexican School of Particles and Fields in October 
2008. 



1 Introduction 

Current evidence points towards the success of the RHIC program at recreating the 
long-sought quark-gluon plasma (QGP) [2], scheduled to be studied in more detail 
at the LHC In the range of accessible temperatures (of order a few times the de- 
confinement temperature), the QGP appears to behave as a strongly-coupled liquid. 
The phenomenon of jet quenching observed at RHIC indicates that partons crossing 
the QGP experience a significant energy loss [3], and so necessitates a framework for 
studying energy loss at strong coupling. 

Motivated by this, a considerable effort has been invested in the study of strongly- 
coupled thermal non-Abelian plasmas by means of the AdS/CFT correspondence 
[UlllEj. At our present stage of knowledge, contact with the real- world QGP produced 
at RHIC or LHC will be feasible only if QCD can be reasonably well approximated 
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by at least one of the various 'QCD-like' gauge theories whose dual description is 
known. In the past few years, encouraging signs in this direction have emerged even 
for the most rudimentary example [I], SU(N) M = 4 super- Yang-Mills (SYM), which 
at zero temperature is completely unlike QCD, but at finite temperature is in various 
respects analogous to deconfined QCD. Studies of energy loss have been conducted 
considering various types of partonic probes of the plasma, including quarks [3 El 
□U], mesons [HI [HI CES 03] , baryons [EH ESI EZj, gluons [H>1 M\ and A;-quarks [IS]. 
For concreteness, in what follows we will restrict attention to the case of a quark. 

The rate of energy loss determined analytically in the seminal works [TJ [HJ, [9] was 
only shown to apply for an isolated quark in the stationary or late-time regimes (i.e., 
when the quark moves at constant velocity, under the influence of an external force 
which precisely cancels the drag force exerted by the plasma, or when the quark has 
decelerated solely under the action of this drag force for a long period of time, and is 
about to come to rest). In the actual experimental setup, however, the configuration is 
neither stationary nor asymptotic: the quark is not externally forced, and slows down 
under the influence of the plasma, whose finite spatial and temporal extent imply that 
the late-time regime is not generically accessible!]] Moreover, the real-world quark is 
not isolated, but is created within the plasma together with an antiquark. In the 
present text we will explore to what extent the rate of energy loss is affected by these 
three issues. 



2 Basic Setup 

For our analysis we need two ingredients: a strongly-coupled thermal SYM plasma, 
and a quark that traverses it. The plasma is known to be described in dual language 
by a Schwarzschild black hole (actually, black threebrane) in asymptotically anti-de 
Sitter space, whose Hawking temperature is identified with T. More explicitly, to 
study Af = 4 SU (N c ) SYM with coupling gyM and temperature T, we must consider 
Type IIB string theory on the (AdS-Schwarzschild)s x S 5 geometry 

R 2 f dz 2 \ 
ds 2 = G mn dx m dx n = — -hdt 2 + dx 2 + — + R 2 dn 5 
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where l s denotes the string length, and Zh is the location of the event horizon. In 
our work the S 5 factor of this geometry will play no direct role. Notice that the field 
theory temperature is identified with the Hawking temperature T# = l/ftZh of the 
black hole. 

It is also known that a static heavy quark corresponds to a static and purely radial 
string extending down from the horizon at z = z^ to a location z = z m where it ends 
on a stack of Nf D7-branes [2D]. These branes cover the four gauge theory directions 
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t, x, and are spread along the radial AdS direction from z = to z = z m , where 
they 'end' (meaning that the S 3 C S 5 that they are wrapped on shrinks down to zero 
size). The D7-brane parameter z m is related to the Lagrangian mass m ^> \/\T of 
the quark through [7] 
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For applications of this formalism to phenomenology, we must choose values of the 
mass parameter z m based on the charm and bottom quark masses, m ~ 1.4,4.8 
GeV. Following [21] . this translates into z m /zh ~ 0.16 — 0.40 for charm and z m /zh ~ 
0.046- 0.11 for bottom. 

The string dynamics follows as usual from the Nambu-Goto action 
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where g a b = d a X m dbX n G mn (X) (a,b = 0,1) denotes the induced metric on the 
worldsheet. 

In the case of vanishing temperature (zh — > oo), we are left in ([T]) with a pure AdS 
geometry, and ([2]) reduces to 

^ / A\ 

Zm = ~ • 4 

2i\m 

One way to verify that the string is indeed dual to a quark is to compute the expec- 
tation value of the gluonic field in the presence of this object. The result is [22] 
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For m — > oo (z m — > 0), this is just the Coulombic field expected (by conformal 
invariance) for a pointlike charge [23]. For finite m the profile is still Coulombic far 
away from the origin but in fact becomes non-singular at the location of the quark. 
We thus learn that the string is actually dual not to a pointlike 'bare' quark, but to 
a 'composite' or 'dressed' quark, including its gluonic cloud. The size of this cloud is 
precisely the length scale z m defined in (T4]), which is therefore understood to be the 
analog of the Compton wavelength for our non-Abelian source. 

Below we will also be interested in the description of a quark-antiquark pair. The 
IIB strings we have introduced above are oriented, and a state with two oppositely 
oriented purely radial strings would correspond to a quark and antiquark that are 
merely superposed. With such boundary conditions, however, the configuration with 
lowest energy is given by a single fl-shaped string that has both of its endpoints at 
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z = z m . This can again be verified by computing the corresponding gluonic profile. 
The result for infinitely massive quarks, at distances \x\ much larger than the dipole 
separation L, is [21] 

CIV^(x) >= 15r »^', (6) 

which has a more rapid falloff than (jHJ), as expected for an overall color- neutral source^] 
An important generic lesson of this type of computations is that the string codifies 
not only the quarks, but also the gluonic field generated by them. Roughly speak- 
ing, the would-be bare quarks/antiquarks are dual to the string endpoints, and the 
surrounding gluonic field is dual to the body of the string. In other words, we have 
learned through AdS/CFT that the usual 'QCD' string exists even for non-confining 
theories, and actually lives in a curved 5 (+ 5)-dimensional geometry. 

Combining the two ingredients above, we know that a quark in a strongly-coupled 
SYM plasma is described by a string on the Schwarzschild-AdS geometry ([!])■ The 
drag force felt by this quark when moving along a trajectory x(t) can then be deduced 
simply by having the string endpoint follow this same trajectory, and determining 
the pull of the body of the string which trails behind it. The works [7J E] derived an 
analytic formula for the corresponding rate of energy loss, 

dE q \ 7T n-„ v 



~lV\T 2 ^L= . (7) 
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As mentioned in the Introduction, this formula was proven to hold only in the sta- 
tionary or late-time regimes, so a priori it is not clear whether it could be relevant 
under conditions analogous to those of the experimental setup. This is the question 
that we now wish to examine. To orient ourselves, we will first inquire about the rate 
of energy loss for an accelerating quark in vacuum, whose dual description involves 
the string moving on pure AdS spacetime. In this case we will have full analytic 
control over the system, which will allow us to develop some intuition on the problem 
at hand. Additionally, our vacuum analysis will provide a useful benchmark against 
which the finite-temperature results can be compared. 



3 Isolated Quark at Zero Temperature 
3.1 Infinite mass 

A quark that accelerates in vacuum would be expected to emit chromoelectromagnetic 
radiation. The first definite characterization of the radiation rate off an accelerating 
quark by means of AdS/CFT was worked out in an important paper by Mikhailov 
|26j . Remarkably, this author was able to solve the highly nonlinear equation of 
motion for a string on AdSs that follows from ([3]), for an arbitrary timelike trajectory 

2 The additional factor of L/\x\ in ((6|) beyond the usual electrostatic dipole behavior L 2 /\x\ e is 
due to the large N c limit [25] . 
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of the string endpoint dual to a heavy quark! In terms of the coordinates used in ([T]) 
(where for now h — 1), his solution is 



dr^(r) 

X^z) = z^p + x»(r) , (8) 

with fi = 0, 1, 2, 3, and x^(t) the worldline of the string endpoint at the AdS boundary— 
or, equivalently, the worldline of the dual, infinitely massive, quark — parametrized 
by the proper time r defined through i] flu x' J 'x l/ = — 1, where x M = dx^/dr. 

Combining ([1]) and ([8]), the induced metric on the worldsheet is found to be 

R 2 R 2 
9tt = -y(z 2 x 2 - 1), g zz = 0, g ZT = -, 

z z z z 

implying in particular that the constant-r lines are null, a fact that plays an important 
role in Mikhailov's construction. In the solution ([H]), the behavior at time t = X°(t, z) 
of the string segment located at radial position z (which according to our discussion 
in the previous section codifies the behavior of the gluonic field at the length scale z 
|27j ) is completely determined by the behavior of the quark/endpoint at a retarded 
time t ret (t,z) obtained by projecting back toward the boundary along the null line 
at fixed t. From the /x = component of (jHD, parametrizing the quark worldline by 
x°(r) instead of r, and using dr = y/1 — v 2 dx°, where v = dx/dx°, this amounts to 

+ tret , (9) 



VT^t 2 

where the endpoint velocity v is meant to be evaluated at t ret . In these same terms, 
the spatial components of (|HJ) can be formulated as 

X(t, Z) = Z V + X (t ret ) = (* - tvct)v + X (t ret ) . (10) 

VI — u 

Working in the static gauge a = t, a 1 = z, the total energy of a string that 
extends all the way down to the boundary — i.e., with z m = 0, corresponding to an 
infinitely massive quark — follows from the Nambu-Goto action as 

-\ 2 
ax ) 



dt I V dz J \ dt J \ dz J \ dt dz J 



Using (Q and fflUj) . Mikhailov was able to reexpress this energy (via a change of 
integration variable z —> t vct ) as a local functional of the quark trajectory, 

y/X f f ,. a 2 — \v x a] 2 



E(t) = — / dt Ict _, ,3 + E q (v(t)) , (12) 
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where of course a = dv/dx°. (A similar expression can be written down for the 
total string momentum.) The second term in the above equation arises from a total 
derivative that was not explicitly written down by Mikhailov, but can easily be worked 
out to be [1] 



which gives the expected Lorentz-invariant dispersion relation for the quark. The 
energy split achieved in (I12p therefore admits a clear and pleasant physical interpre- 
tation: E q (associated only with information of the string endpoint) is the intrinsic 
energy of the quark at time t, and the integral over t ret (associated with the body 
of the string) encodes the accumulated energy lost by the quark to its gluonic field 
over all times prior to t. No less remarkable is the fact that the rate of energy loss 
for the quark in this strongly-coupled non-Abelian theory is found to be in precise 
agreement with the standard Lienard formula from classical electrodynamics! The 
AdS/CFT correspondence thus teaches us that, in this very unfamiliar non-linear 
setting, the energy loss turns out to depend only locally on the quark worldline. This 
feature has been argued in [28J to lead to an upper bound for the energy of a quark 
at finite temperature. 

3.2 Finite mass 

It is interesting to consider how these results are modified in the case z m > 0, where 
the mass m of the quark given by (j4j) is large but not infinite. In this case, we ought 
to impose boundary conditions on the string not at z = but at z = z m : given 
the worldline x(t) of the quark, we must require that the string worldsheet satisfy 
X(t, z m ) = x(t). Moreover, we need only determine the behavior of the string in the 
region z > z m . 

Our task is to reexpress (jSJ) in terms of the data x(t) at the new boundary z = z m , 
instead of the (now merely auxiliary) data at the AdS boundary 2 = 0, which we 
will henceforth distinguish with a tilde: x(t). For simplicity, we will carry out this 
translation explicitly only in a setup where the quark moves purely along direction 
x = x 1 , which is all that we will need for our analysis in subsequent sectionsJl 

As was explained in [T], it is possible to rewrite the velocity v and acceleration a 
of the 'auxiliary endpoint' at z = in terms of the velocity v = dx/dt = dtX(t,z m ) 
of the actual string endpoint and the momentum density IPj , which, when evaluated 
at z — z m , controls the external force F = (^/A/27^)^*(i, £ m ) acting on the string 

3 The corresponding formulas for arbitrary motion can be found in [29j and a more detailed 
explanation in |30j . 
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endpoint, or equivalently, on the quark. One finds 



i — z^v n 
n (i-, 2 )^(i-4n 2 )^ 
(1 - 
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where we have abbreviated II = ITj. Using f|T2|) (decorated with appropriate tildes), 
we can then express the total energy of the string at time t in the form 
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As before, the first term represents the accumulated energy lost by the quark at all 
times prior to t: it is the generalization to the m < oo case of the Lienard formula 
f|T2|) deduced by Mikhailov (indeed, for m — > oo the integrand in (I15p reduces to the 
Lienard expression for motion in one dimension, (\^\/2tt)F 2 /m 2 = (v / A/27r)a 2 /(l - 
t> 2 ) 3 ). The second term again denotes a surface term and gives the modified dispersion 
relation for the finite-mass quark, 
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Similar expressions can be written down for the momentum [T]. 

Notice that ( {TBI) differs from the result expected for a pointlike quark, E = 7m. 
This reflects the fact that the fundamental source dual to a string that terminates 
at z m > is indeed not pointlike. According to the standard dictionary mapping the 
radial coordinate in the curved geometry to a length scale in the gauge theory |27j . 
it has a linear size of order z m . It is only because of this extended nature that, as we 
saw above, to characterize its state one needs to specify not only the velocity v but 
also the parameter F (or d z X(t, z m )) that encodes its shape. The crucial point here 
is that, as we already noted below ([5]), the source in question should not be thought 
of as a bare quark, but as a 'dressed' or 'constituent' quark, surrounded by a gluonic 
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cloud with thickness z m [ 

Another salient feature of the energy of the quark given by the expression 
the fact that it diverges as the value of the external force approaches 



is 



crit 



(17) 



The reason for this is easy to understand on the string theory side. To exert a force 
F on the string endpoint, within the D7-branes we must turn on an electric field that 
has strength F m = F at z = z m . The physical origin of the bound on the force is 
the fact that, for F i > Fq"* = y/X/2irz^ n , the creation of open strings is energetically 
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favored, and so the system is unstable. According to (fl6l) . then, the energy of the 
constituent quark diverges precisely at the point where the external force becomes 
capable of nucleating quark-antiquark pairs. 

To summarize, from our study of a quark accelerating in vacuum we learn that: 
i) the total energy of the string includes not only the intrinsic quark energy, but also 
the accumulated energy previously radiated by the quark; ii) the quark dispersion 
relation arises from a surface term on the string worldsheet; in) the radiated energy 
depends localy on the quark worldline; iv) when the quark has a finite mass, it is 
automatically non-pointlike, and there are consequently non-trivial modifications to 
its dispersion relation and radiation rate. 

Before closing this section, we would like to make an additional observation. It is 
natural to expect the energy split achieved in (1121) or (1151) to be somehow reflected 
in the geometry of the string worldsheet. It is interesting then that, on the string 
worldsheet dual to a quark undergoing arbitrary accelerated motion, there appears a 
black hole [I], whose event horizon naturally divides the worldsheet into two causally 
disconnected regions (the appearance of a worldsheet black hole had also been noted 
previously at finite temperature [32], [33]). It seems plausible then to interpret the 
regions outside and inside the black hole as corresponding respectively to the quark 
and the gluonic field, and postulate that the energy flow across this horizon should be 
related to the energy radiated by the quark [Mill]. This possibility has been studied 
more closely in [35], [361 EZ] • 

4 Isolated Quark at Finite Temperature 
4.1 Constant velocity 

Having understood the rate of energy loss for a heavy quark that moves in the SYM 
vacuum, in this section we restore < oo — and consequently h < 1 — in the 
metric (0Q), to study the same quantity in the case where the quark moves through 
a thermal plasma. A thorough generalization of Mikhailov's analytic results [26] to 
this finite temperature setup would require finding the exact solution to the Nambu- 
Goto equation of motion for the string on the AdS-Schwarzschild background, for 
any given trajectory of the string endpoint at z m > 0. Sadly, we have not been 
able to accomplish this feat. Nevertheless, based on the results discussed in the 
previous subsection, we expect the total energy of the string at any given time to 
again decompose into a surface term that encodes the intrinsic energy of the quark 
and an integrated local term that reflects the energy lost by the quark. 

There are two easy cases where one can show analytically that this expectation is 
borne out. The first is the case of a static quark, where there is of course no energy 
loss, and the string energy is correspondingly given by a pure surface contribution, 
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which we naturally interpret as the rest mass of the quark in the thermal medium. 
An important difference with respect to the T = case analyzed in the previous 
subsection is that here the surface contribution arises not only from the lower (z = z m ) 
but also from the upper (z = Zh) endpoint of the string. This is in fact the generic 
situation in the T > (zh < oo) case. Since Zh marks the position of an event 
horizon, for any finite coordinate time t the value of the surface contribution at Zh is 
not influenced by the behavior of the z < Zh portion of the string, but depends only 
on the string's configuration at t — > — oo. The same interpretation can be then carried 
over to the gauge theory: for arbitrary finite-temperature configurations, the terms 
in the quark dispersion relation arising from the upper string endpoint will encode a 
contribution to the energy of the state that depends solely on the initial configuration 
of the quark+plasma system. Throughout the evolution, causality guarantees that 
the behavior of the SYM fields at spatial infinity can only be affected by the initial 
configuration at t — >• — oo, so we can equivalently think of the surface terms at Zh as 
encoding information on the asymptotic boundary conditions for the system. Indeed, 
for dynamical processes, the radial location z — Zh'va. AdS-Schwarzschild corresponds 
to the deep IR of the gauge theory. 

The only other finite-temperature solution that is known analytically corresponds 
to a quark moving at constant velocity v, and takes the form [7J [8] 



X(t,z) 



' z h ( z h + z\ z h _ x ( z 
t -In H tan — 



Zh-z) 2 \z h 



(19) 



Imitating the procedure of the previous section, we can rewrite this solution by slicing 
it along null worldsheet curves, 

X ^ z ) = (i - ,,2)1/4 tan ~ [ Zh (i- v 2y/4 ) + > ( 2 °) 

allowing the total energy of the string to be split into [1] 




As expected, the integrated term in the top line recovers the result for the stationary 
rate of energy loss (171). The terms in the second and third line, then, codify the energy 
Eg intrinsic to the quark, including, as we discussed above, contributions dependent 
on the initial or boundary conditions. The total momentum of the string can be 
processed similarly. 
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4.2 Accelerated motion 



Having gained some intuition from the analysis of a quark moving at constant speed 
relative to the strongly-coupled plasma, let us now turn our attention to the more 
general situation where the quark accelerates. For concreteness, we will consider the 
early-time behavior of a quark initially at rest, which is accelerated along a line by 
a time-dependent external force that is turned off after a short interval, allowing the 
quark to move thereafter only under the influence of the plasma. 

In [1] we integrated numerically the Nambu-Goto equation of motion for the dual 
string. Our results show a qualitative difference between the initial stage (0 < t < 
Release) where the quark is accelerated by means of the external force F(t), and the 
second stage (t re iease < t < t ma , x ) where it moves only under the influence of the plasma. 
In the former stage, the rate of energy loss, for values of m in the neighborhood of 
the charm or bottom masses, is in fact nearly identical to the corresponding rate in 
vacuum, given by the modified Lienard formula (fT5l) . In what follows we will elaborate 
a bit on the behavior in the latter stage, which would appear to be more relevant 
from the phenomenological perspective. A more detailed discussion of the energy lost 
by the quark in both stages can be found in [T]. 




Figure 1: Quark velocity as a function of time (in units of 1/vrT) from our numerical 
integration (in red) compared against ( |25i) with the value of /i deduced in [7] (in 
black), and (I25I) with \x chosen to fit the data (in light blue), for a) z m /zh = 0.2 and 
b) z m /z h = 0.4. 



The most direct way to inquire whether the output of our numerical integration for 
the accelerated quark conforms to the the constant-velocity (or late-time) results of 
[7J [S] is to compare the corresponding quark trajectories. If one assumes a dispersion 
relation p q = jM^v as in [7J (which is what the vacuum expression (fT6l) predicts in 
the unforced case), then the rate of energy loss ([7]) is equivalent to the statement that 
the friction coefficient 

^-1%. (22) 
p q at 
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takes the constant value \i = ny/XT 2 /2M^ ia . The corresponding equation of motion 
for the quark is then 

^ = -Hl-« 2 ), (23) 



whose solution is [7J 



v(t) = Vvelease — . (24) 

3 2/i(t-i release ) 



^Areleasc + U ^release)' 



In our numerical results f re iease <C 1, so we are only able to test the non-relativistic 
version of (1241) . 

v(t) = w"" (t "^« ) (25) 

(and, since our integration is limited to small time intervals, we would in effect see just 
the linear portion of this function). A comparison between this analytic prediction 
and our numerical results for t > t TC \ casc is given in Fig. [TJ It is evident from this plot 
that, in the early stage of motion covered by our analysis, the quark dissipates energy 
at a rate much lower than the late-time result of [7J. Indeed, for z m /zh = 0.2, 0.3, 0.4 
the late-time friction coefficient is respectively fii atc /iiT = 0.25,0.41,0.59, but our 
numerical results for v(t) are best approximated by /Wiy/ 71 "^ = 0.08,0.15,0.26. 

The relation between energy loss and the appearance of a black hole on the world- 
sheet in this finite temperature setting had been noted already in [321 E3] for the case 
of a stationary configuration. The non- stationary case and the parallel with the zero 
temperature situation were discussed in [Tj |35| [36]. 

An interesting prediction of AdS/CFT is the existence of a subluminal limiting 
velocity v m = y/l — z^/z^ ~ 1 — (\^XT/2m) 4 for the quark traversing the plasma. 
This follows simply from the fact that the quark velocity is dual to the coordinate 
velocity of the string endpoint, and v = v m corresponds to a proper velocity at z = z m 
equal to that of light [38] • A more general bound involving the external force F can 
be derived by requiring the Nambu-Goto square root to remain real p], which reduces 
to v < v m when F — 0. This same restriction on the velocity can be seen to arise 
from microscopic calculations of the meson spectrum [39l 00] . Its validity for isolated 
quarks has been emphasized in [4~T] . 



5 Pair Creation within the Plasma 

As was mentioned in the Introduction, in the experimental setup the heavy quark is 
created within the plasma together with its corresponding antiquark, and the presence 
of the latter would be expected to substantially modify the gluonic fields in the vicinity 
of the quark, consequently affecting its evolution. In [I], building upon the previous 
work [7J, we have numerically studied the evolution of a heavy quark and antiquark 
that are created within the plasma at time t — and then separate back to back. 
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In dual language, this corresponds to a string with both of its endpoints on the D7- 
branes at z = z m , such that the endpoints are initially at the same spatial location 
but have initial velocities in opposite directions. 

In setting up the problem, one realizes that on the gauge theory side there are 
actually two distinct quark-antiquark configurations: the product of a fundamental 
q and an antifundamental q can lead to a q-q pair either in the singlet or the adjoint 
representation of the SU(N C ) gauge group. Interestingly, we find that there is a 
counterpart to this in the gravity side, because there are precisely two distinct types 
of consistent initial conditions for the string. The most obvious possibility is a string 
that is completely pointlike at t = j7], and then expands into a fl-shape as its 
endpoints separate. Through the standard recipe for correlation functions [5], it is 
clear that, initially, this configuration sets up no long-range chromoelectromagnetic 
field, and so describes the singlet. The only other allowed possibility is a string that 
is initially linelike [I], extending from z m all the way up to Zh, and expands for t > 
into a A-shape. Such linelike string is precisely the system considered in (121 [15] to 
model a color source in the adjoint representation. Due to its extended nature, it sets 
up long-range supergravity fields that translate into a long-range gluonic field profile, 
indicative of a source with net color charge. Curiously, we find that the initial quark 
velocity is freely adjustable in the adjoint, but not the singlet, configuration — in the 
latter case it is invariably fixed at the limiting value v = v m . 

As the quark and antiquark separate, we expect them to eventually enter the 
late-time regime where their rate of energy loss is given by ([7]) and their trajectory 
is given by ( 1241) (with Release and v release regarded as adjustable parameters). The 
agreement between the analytic and late-time numeric results is most cleanly seen 
if instead of comparing graphs of x(t) or v(t) for the quark (where one would need 
to look at t — > oo), one examines the plots of v(x), where the analytic behavior for 
constant /i takes the simple form 



which is linear with slope —\i near the final rest point x = x^ (whose value is meant 
to be adjusted to fit the data). These plots are shown in Fig. [21 for mass parameter 
Zm/zh = 0.2 (in the neighborhood of the charm quark), corresponding to a limiting 
velocity v m = 0.9992. It is evident from the figure that the late-time behavior is 
well- described by (126!) in all cases. 

In Fig. [21 we see that there is an initial period where the behavior differs from 
the late-time frictional evolution (I26p . This difference is clearly more significant for 
the singlet than the adjoint case. For values of the mass in the neighborhood of the 
charm or bottom masses, our numerical results indicate that the initial evolution is 
essentially identical to what it would be in the absence of the plasma. This picture 
seems rather close to the phenomenological discussion given in [19j (in the context 
of collisional energy loss): when the singlet quark-antiquark pair is formed within 
the plasma, there is a delay before the interaction between the newly created sources 
and the plasma can set up the long range gluonic field profile that is responsible for 




(26) 
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Figure 2: Quark evolution (velocity as a function of traveled distance, in units of 
1/irT) for five different initial conditions. The dotted curves show the results of 
our numerical integration, contrasted against fits in solid red that use the analytic 
expression (126]) . with the value \x = 0.25 obtained in [7] and an optimal choice of the 
stopping distance x^. The three dotted curves starting at v = v m describe singlet 
configurations with different initial velocity profiles and energies. Notice in particular 
that the purple curve describes a situation where the quark and antiquarks turn 
around and come to rest while approaching one another. The two remaining curves 
arise from adjoint configurations with different energies and initial quark velocities. 
The vertical arrows mark the transition distance where each trajectory enters the 
late-time regime described by (I26p . 



the late-time dissipation. To examine in more detail the transition to the late-time 
behavior (|26|) . for a variety of trajectories we can determine the point (x/, Vf) beyond 
which the numeric v(x) curve agrees with the analytic curve (|26|) to a given accuracy 
/. (In pQ we considered / = 5 or 10%.) The results are schematically indicated by 
the arrows in Fig. [21 

By the time when the quark moves beyond Xf{vf) and therefore enters the late- 
time regime, it is certainly insensitive to the presence of the antiquark. It is natural 
then to ask how the transition distance Xf(vf) compares against (half of) the length 
L max (t)) beyond which the quark and antiquark are screened from each other by the 
plasma. Based on what we just said, we know that for a given velocity we must have 
-^max/2 < Xf, but the actual comparison will tell us whether the transition to the 
regime where the quark experiences a constant drag coefficient occurs right after the 
quark and antiquark are screened from each other by the plasma, or if there is an 
intermediate regime where the quark moves independently from the antiquark but 
nevertheless feels a drag force that differs from the stationary result of [7J [H [9] , as we 
found when applying an external force to the isolated quark in Section 14.21 

The screening length for infinitely massive quarks in M = 4 SYM was computed 
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in [T2l [T3] . by considering a quark- ant iquark pair moving jointly through the plasma 
(a related length was plotted in [H]). Over the entire range < v < 1 its behavior 
may be approximated as 



0.865 
ttT 



(27) 



while in the ultra- relativistic limit, it can be shown analytically that 



^3-3/44^/2 
vrT T(l/4) 2 



0.743 



1 - v 2 ) l ' A for v -»■ 1. 



(28) 



The full curve L max (t>) does not deviate far from this asymptotic form, so a decent 
approximation to it is obtained by replacing 0.743 — > 0.865 in ( 1291) . to reproduce the 
correct value at v = (at the expense of introducing a 16% error as v — > 1) [T2] : 



0.865 

iff' 



,2x1/4 



(29) 



A comparison between the two approximations ( 1271) and ( 1291) is shown in Fig. [3j 
overall, the exponent 1/3 is better than 1/4 in the sense that it implies a smaller 
squared deviation from the numerical results, even though 1/4 leads to a smaller 
percentage error in the range v > 0.991 (7 > 7.3). An attempt to better parametrize 
the deviation away from the ultra-relativistic behavior was made in [13]. In any case, 
one should bear in mind that the region of principal interest at RHIC is not really 
v —> 1, but jv ~ 1. 




Figure 3: Screening length L max (in units of l/2irT) as a function of velocity (in 
black) compared against the approximations ( 1271) (in red) and ( 1291) (in blue). 



To compare against the transition distance Xf defined above, in pQ we extended 
this calculation to the case of finite mass. Fig. H] shows the result for a mass in the 
neighborhood of the charm value. The main difference with respect to the m — ► 00 
case is that now the 'ultra-relativistic' region would refer to the limit where the pair 
velocity approaches the limiting velocity v m < 1. The screening length L max (t>) is still 
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relatively well approximated in the full range < v < v m by the natural modification 
of the z m = fit (1271), 

W*)*-^^-* 3 ) 173 . (30) 



ttTv 



This approximation becomes worse as z m /zh is further increased. (In all cases, the 
fit analogous to (1291) . 

(31) 



0.865 



2\l/4 



does a poorer job than ( 1301) .) The behavior in the v — > u m region can still be deter- 
mined analytically, and turns out to be 



v 2 — v 2 

m 



1 zl[l - {z m / z h f - v 2 



max 2nT v m (l-v 2 m y/* 2vrT 4^1 _ (2m /^)4 ' 
where we see that the 1/4 exponent in (12"B|) changes to 1 at finite quark mass. 



(32) 





0.992 0.99 



Figure 4: (a) Screening length L max (in units of 2irT) as a function of velocity for 
z ml Zh — 0-2 (in black) compared against the z m > fits (130|) (in red) and (13T|) (in light 
blue), (b) Expanded version of the same plot, showing that (130|) gives a relatively 
good approximation up to velocities very close to v m , where the asymptotic linear 
behavior (13"2"1) (dotted dark blue) sets in. 



We can now finally perform the promised comparison between the transition dis- 
tance x/(v) and (half of) the screening length L max (v). The result, for singlet config- 
urations, is shown in Fig. [51 where we see that the two separations are of comparable 
magnitude and scale with velocity in a similar mannero A similar agreement was 
reported in [45] (later related work may be found in [16]). Notice that this is in spite 
of the fact that the two relevant string configurations are quite different, with motion 
respectively along and perpendicular to the plane in which the string extends. 

We conclude then that the transition to the constant-drag-coefficient regime takes 
place immediately after the quark and antiquark lose contact with one another. That 



4 For adjoint configurations, the two lengths also agree but have a trivial v depencence: Xf is seen 
in Fig. [2] to be essentially equal to zero, and the adjoint quark-antiquark potential is suppressed by 
a factor of l/N% (see, e.g., [44]), implying that L max = at large N c . 
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Figure 5: Transition distance Xf(v), with z m /zh = 0.2, for / = 0.1 (red) and / = 0.05 
(magenta), compared against (half of) the screening length L max (v) (green). The 
vertical axis is in units of 1/ttT. 



is to say, unlike what we found for the forced isolated quark in Section 14. 2\ here there 
is no intermediate stage where the quark and antiquark decelerate independently from 
one another at a rate that differs substantially from the late-time result of [3, [8]. 

To summarize, in the more realistic scenario where the quark and antiquark are 
created within the plasma and then separate back-to-back, we find that there are two 
distinct stages of the evolution. At the beginning, the quark essentially does not feel 
the plasma and slows down only under the influence of the antiquark, just as it would 
in vacuum. To the extent that this result, with all its simplifying assumptions, might 
conceivably be extrapolated to the experimental context, this initial stage would not 
differ significantly between heavy ion and proton-proton collisions. Later, there comes 
a point when the q-q separation L qq - becomes large enough, and the quark velocity 
low enough, that the pair crosses the relevant screening length, i.e., L qg - « L max (v). 
Beyond this point the quark starts feeling the plasma but no longer feels the antiquark, 
and evolves exactly as the expressions (ED) and (124j) predict. This then extends and at 
the same time delimits the region where the analytic results of [3, El E] can justifiably 
be used to model energy loss in heavy ion collisions. 
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